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We give an exposition and numerical studies of upper hedging prices in multi- 
nomial models from the viewpoint of linear programming and the game-theoretic 
probability of Shafer and Vovk. We also show that, as the number of rounds goes 
to infinity, the upper hedging price of a European option converges to the solution 

of the Black-Scholes-Barenblatt equation. 
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O 1 Introduction 

O' 

The Black-Scholes formula for the geometric Brownian motion model and the Cox-Ross- 
Rubinstein formula for the binomial model are now treated in many standard textbooks 
/\^ • on mathematical finance (e.g. [131 1211 [22] )• Since these models are complete, the exact 

j^ ■ price for any contingent claim is determined by arbitrage argument. On the other hand, 

incomplete models such as the trinomial model are only briefly mentioned in the textbooks 
because of difficulty associated with indeterminacy of prices of contingent claims. 

In fact only a few explicit results seem to be known on upper hedging prices for the 
discrete-time multinomial models. The purpose of this paper is to give an exposition 
of the exact and the asymptotic behavior of upper hedging prices of contingent claims 
in multinomial models. We also show that the asymptotic upper hedging price of a 
European option is described by the Black-Scholes-Barenblatt equation (e.g. [1], [26], [9], 
[12] . Chapter 4 of [IS]). The Black-Scholes-Barenblatt equation is usually considered for 
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uncertain volatility models in continuous time and its relation to multinomial models does 
not seem to be stated in literature. 

The advantage of discrete-time multinomial models is that we can exactly compute the 
upper hedging price by backward induction for moderate number of rounds and various 
approximations to the upper hedging price can be compared to the exact value. Basic 
facts on the upper hedging price for discrete models are well explained in Chapter 4 of the 
first edition of Musiela and Rutkowski [13]. General treatments of incomplete markets 
are given in Chapter 5 of [TO] and [18]. However they are concerned with continuous 
stochastic processes and do not contain much numerical information on the behavior of 
upper hedging prices for discrete models. 

An extensive numerical study of hedging in incomplete markets is given in [3]. Its 
authors consider a hedging strategy which minimizes the mean-squared error to the payoff 
of a contingent claim under Markov-state dynamics. As we see in Section [2] below, for 
studying the behavior of upper hedging prices we can not make convenient stochastic 
assumptions such as the Markov property. Results more relevant to the present paper 
have been given in [T7] and [S] by convex ordering argument. In particular for discrete 
time models with bounded martingale differences, [17] proves that the upper hedging price 
of a convex contingent claim is given by the extremal binomial model. We reproduce this 
fact in Proposition 12.41 below by linear programming argument. 

In this paper we use the framework of game-theoretic probability by Shafer and Vovk. 
We prefer the framework because of the following reasons. First, in game-theoretic proba- 
bility only the protocol of a game between "Investor" and "Market" is formulated without 
specification of a probability measure. Therefore there is no need to distinguish a risk 
neutral measure from an actual (or a physical) measure and to consider the equivalence 
between them. This is advantageous because the extremal risk neutral measure corre- 
sponding to the upper hedging price usually has a support smaller than those in the 
interior of the set of risk neutral measures. Second, some strong properties of a price path 
of Market can be proved in game-theoretic probability without any stochastic assumption 
(e.g. see [25], [27], [28] and references therein). As a referee pointed out [^ studies non- 
probabilistic approach for pricing in continuous time. Third, the notion of upper hedging 
price is of central importance to game-theoretic probability as shown in our recent works 
( |24] . [20] ) on game-theoretic probability. 

The organization of the paper is as follows. In Section [2] we give a linear programming 
formulation of upper hedging prices in multinomial models and state some basic facts in 
several propositions. In Section [3] we give some simple bounds for upper hedging prices. 
Then in Section Hj we show that, as the number of rounds goes to infinity, the upper 
hedging price of a European option converges to the solution of an additive form of the 
Black-Scholes-Barenblatt equation. In Section |5] we present numerical studies on the 
accuracy of the partial differential equation and other approximations. Some concluding 
remarks are given in Section |6l 



2 Formulation of upper hedging price 

In this section we formulate the upper hedging price for a multinomial game from the 
viewpoint of linear programming and the game-theoretic probability of Shafer and Vovk 
[19] . Also we show some known facts on upper hedging prices. 

Let Af C M be a finite set containing both negative and positive elements. The protocol 
of the multinomial game of N rounds with the initial capital oUCq = a is written as follows. 

K.Q = a 

FORn= 1,2,..., AT. 

Investor announces M„ G M. 

Market announces x^ G A". 

/C„ := )Cn-l + MnXn- 

END FOR 

We call X'^ the sample space and ,^ = xi . . . xn G '^^ a path of Market's moves. For 
1 < n < N , ^^ = Xi . . . Xn & X"- is a partial path. Investor's strategy P is a function 
specifying M„ based on ^""^ 



Xl ■ ■ ■ Xji—l- 

V -.xi... 



x„_i H-). M„. 



with some initial value Mi = 'P(n), where D denotes the initial empty path. When 
Investor adopts V, his capital at the end of round n is written as a + /C^, where 



K = Y.'^{r')x^ 



We can write the progression of the game in a rooted tree with the root D. For n < N, 
^^ is an intermediate node branching to ^"x„+i, x„+i G X. The final nodes are ^ G X'^ . 

We call a function / : X^ — )■ M a payoff function or a contingent claim. The upper 
hedging price (or simply the upper price) of / is defined as 

Eif) = inf{a I 3V s.t. a + /C^(0 > /(O, V^ G X'^} 

and the lower hedging price is defined as 

E{f) = -E{~f). (1) 

The upper hedging price and the lower hedging price are often called seller's price and 
buyer's price, respectively. A strategy V with the initial capital a satisfying a + /C^(^) > 
/(Oj^^ ^ ^'^y is called a superreplicating strategy for /. 

The problem of obtaining the upper hedging price can be formulated in linear pro- 
gramming. Let X = {oi, . . . , Ofc}. For the single step case A^ = 1, E{f) is obtained as the 
following minimum: 
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For the two-step case N = 2, Investor can choose his investment at round 2 depending 
on the Market's move in the first round. Therefore E{f) is written as the following 
minimum: 
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For general A^, the coefficient matrix A^ : A;^ x (1 + {k^ — l)/{k — 1)) is recursively 
defined as 

/ /ai\ \ 
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where 1^ is a fc-dimensional vector of I's, A^v is A^v without the first column and 
denotes the Kronecker product. Then 

E{f) = min(l ... 0) m, s.t. Antti > /, 

where / is a fc^-dimensional vector consisting of f{^),^ € X^, and 



(3) 



m^ = (a Ml Maiai . . . M 
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Since the size of An grows exponentially with N, it becomes difficult to directly 
solve the linear programming problem for a general path-dependent contingent claim /. 
Exploiting the recursive structure of the coefficient matrix A^, the linear programming 
problem can be solved by backward induction, i.e. by solving the single-step optimizations 
for N,N — !,...,!. This will be explicitly described in ([H]) below. Therefore the single- 
step optimization in ([2]) is essential. However even with backward induction, for a general 
path-dependent /, the number of single-step optimizations grows exponentially with N. 
This is because the number of nodes of the game tree grows exponentially with A^ and 



the single-step optimization for the backward induction is performed at each node of the 
game tree. This difficuhy is somewhat mitigated in the case of a European option /, 
where / is a function oi S^ = Xi + ■ ■ ■ + x^ only. Then the game tree can be collapsed 
according to the values oi Sn = Xi + ■ ■ ■ + Xn and the number of nodes of the collapsed 
tree grows only polynomially with A^. In fact, for generic values of ai, . . . , a^, the number 
of values taken by Sn is ( ^_i )■ 

Now we consider the single-step optimization in ([2]). Let X contain / negative elements 
and m = k — I positive elements, which are ordered as 

> ai > a2 > ■ ■ ■ > a^ , < a^^ < a^ < ■ ■ ■ < a^. 

Note that we allow the case a^ = 0, which needs some special consideration. For the 
single-step game the following result is given in Proposition 4.1.1 of [13]. We give our own 
proof based on consideration of dual linear programming problem to ([2]). 

Proposition 2.1. The upper hedging price of f in the single-step game N = 1 is given 
by 
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and consider the following dual problem to ([2]): 
• • /(^m)) P -^ max 

Note that the coefficient matrix is 2 x /c. Therefore the maximum is attained by a basic 
solution involving two variables. 

First consider the case a^ > 0. If we choose two variables either from {pi ,P2 , ■ ■ ■ ^PT} 
or from {pf,p2,.-.,p^}, then the solution does not satisfy p > and it is infeasi- 
ble. Therefore for a feasible basic solution we need to choose one variable p^ from 
{Pi,P2, ■ ■ ■ ,PT} s-^d another variable p^ from {pt,pty ■ ■ ■ ^Pm}- Then the basic solution 
is given by 

Pi = Pr.j = + Pj = Pj;i = + ' + (6) 

and the value of the objective function is given by {a^ f{aj) — a^ f {a^)) / {a^ — aj). 

Now consider the case a^ = 0. Then for j = 1 we have {a^ f{a^) — aJ f{a^))/{a^ — 
a~) = /(O), which is consistent with (jS]) and ([H]). Therefore this case does not require a 
separate statement. This proves the proposition. D 

By ([T]) the lower hedging price E_{f) for the single-step game is given as 

E{f) = min ' ^'^ ' ' '•"• ^ 



hj \ aj — a^ 



The vector p satisfying ([5]) is a probability vector over X such that its expectation is 
zero: 

/ m 

i=l j=l 

Let Vi denote the set of probabihty vectors satisfying (j^]) . p G Pi is called a risk neutral 
measure on X. From (^ we have E{f) = maxpg-p-^ E^lf], where -E^[-] denotes the usual 
expectation under p. This duality is well known in great generality (see Chapter 5 of 
[TOj ) . However the explicit expression of the upper hedging price in (^ is useful for the 
purpose of numerical backward induction for N >2. 

Remark 2.2. We can also consider the dual problem to ^ for the N-step problem. It is 
easy to see that the set Vn of risk neutral measures p = {p{^),C, G X"-} is characterized 
as follows: 

PEVn <^ p{xi...xn) =Pi{xi)y~p2{x2\xi)y~ ■■■ y~pN{xN\^^~^), (7) 

s.t. p,{- 1 e-1) e Pi, J = 1, . . . ,n, ve'-' G x^-\ 

Also by duality we have E{f) = maxpg-p^ -S^[/]- Note that ([7]) holds even if p{C,) = for 
some ^ and some conditional probabilities on the right-hand side are not defined. At each 
node of the game tree, the conditional distribution of the extremal risk neutral measure 
corresponding to E{f) is given by a basic solution of the form (jlj). 

Based on (jl]), the backward induction for obtaining E{f) for the A^-step game is 
described as follows. Define /(-, N - n) : X"" ^ R, n = N, N - 1, . . . ,0, by 

'atf{Ca^,N-n-l)-arf{Cai,N-n 
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with the initial condition /(^^, 0) = /(^^), C,'^ G X'^ . Summarizing the arguments above 
we have the following proposition. 

Proposition 2.3. The upper hedging price of f in the N-step game is given by 

E{f) = f{n,N). 

Consider the case that the payoff / is defined on the whole M^. Then we have the 
following result. 

Proposition 2.4. Suppose that f is convex. Then E{f) is given by the binomial model 
concentrated on two outermost values {a^,a^}. For the case af>0 E_{f) is given by 
the binomial model concentrated on two innermost values {a]",a^}. For the case a^" = 
^(/) = /(0). 

This result was proved by Riischendorf ( [17] ) by a convex ordering argument in a more 
general setting, where A" is a bounded interval. See also [5]. Corresponding results for 
concave / are obtained by using ([T]). We now give a simple proof of this proposition using 



Proof. First consider the one-step game A^ = 1 . By convexity of / it is easy to check that 
for any a^ and a^ the following inequalities hold: 

<fM-arf«) ^ 4fioi)--^fi4) ^ p^^("g::^^("^) ifa|>0, 
am-a« 4~"i 1/(0) ifa^ = 0. 

Hence the proposition holds for the single-step game. 

For A^ > 1 we can use the induction. We note that f{C,^, N—n) is a convex combination 
of /(^"x„+i . . . xjy), Xn+i . . .x^ E X^^"", where the weights of the combination are given 
by the binomial model and they do not depend on C,"'- In particular /(^", N — n) is convex 
in Xn- □ 

Remark 2.5. As pointed out in Remark 2 of [F?^ , the above result holds if f is component- 
wise convex in every Xi, i = 1, . . . ,N. 

3 Some bounds for upper hedging prices 

In this section we present some simple inequalities for upper hedging prices. 

We first investigate the relation between binomial model and general multinomial 
model. Pick a negative element a^" and a positive element a^ from X and restrict the 
move space X of Market to the two element set Xij = {a~,a^}. Then the multinomial 
game is reduced to a binomial model, where the price of any contingent claim is given by 
an arbitrage argument. Let Eij{f) denote the price of / under the binomial model A^j. 
For the two-element set Xij there is no need to take the maximum in (IHl). On the other 
hand, for a multinomial model X we take the maximum in ([8]) at each node of the game. 
This gives the following lower bound for E{f): 

maxE,j{f)<E{f). (9) 

Similarly for the lower hedging price we have 

minE,,,(/)>E(/). 

The above inequalities can be generalized by considering nested move spaces of Market. 
Explicitly writing X in the multinomial game, let E{f \ X) and E_[f \ X) denote the upper 
hedging price and the lower hedging price of multinomial game with X as Market's move 
space. Consider two nested move spaces X G X' oi Market. Then the same consideration 
as above gives the following inequalities: 

m I X') < E{f I X) < E{f I X) < E{f I X'). (10) 

In Section [5] we compare upper and lower hedging prices in trinomial and quadnomial 
models. 



We can also consider dynamic restrictions of the move space of Market. For example, 
we can consider the maximization in ([8]) only for even n = 2h and use the maximizer i*,j* 
from this round for the subsequent round n = 2h + 1. Or we can maximize e.g. every 
10th step. By this dynamic restriction we again have a lower bound for the upper hedging 
price. As we discuss in the next section, this pruning of maximizations is conceptually 
close to discretization of partial differential equation in ( !T3|l . 

We now present another bound when / is a European option depending only on Sn- 
Let / be defined on the whole M. Assume that / has the first derivative /' which is of 
bounded variation. Then /' is written as a difference of two non- decreasing functions 
(Section 5.2 of J16j). By taking the indefinite integral of /' we see that / is written as a 
sum of a convex function and a concave function: 

/ = /i + /2, /i : convex, /2 : concave. 

By the subadditivity of the upper hedging price (Section 8.3 of [IS]), we have 

E{f)<E{f,) + E{f,). (11) 

The bounds in iQ and flTTl) are very simple. Unfortunately as seen from numerical 
examples in Section [5] these bounds are generally not very tight. 

4 Limiting behavior of upper hedging price of an Eu- 
ropean option 

In this section we derive the limit of the upper hedging price of a European option as 
A^ — )■ cxD in an appropriate sequence of games. Let C^ denote the space of functions 
M — )• R with compact support and continuous second derivatives. Let f E C^. We 
consider a sequence of multinomial games with A^ rounds, where the payoff f^ for the 
A^-th game is given as 

fN{^'') = f{SN/VN), SM = x, + --- + xr,. (12) 

Note that the expected value under the two-point distribution in ([6]) is zero and the 
variance is given as 
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In view of this define the maximum variance and the minimum variance of Vi as 

-2 - + ^ 2 - + 

a = —a^ a^i > g_ = —a^ a{. 
We now state the following theorem. 



Theorem 4.1. Let f E C^ and let f^ be defined by [W) . Assume q^ > 0. Then 

lim ^(/Ar)= 0(0,1), 

where (j){s,t), s G M, < t < 1, satisfies the following partial differential equation 

^A.( f^ ^"^"m ,^ j^' = ^^ ifft0(.,t)>O, 

7nns,t) = —j—(p{s,t), < .r% 1, ,x n ^3 

with the boundary condition (j){s, 0) = /(s), s G M. 

Similarly the following partial differential equation describes the limiting lower price 
of /at. 

^^( ,^ ^'^\( ,^ j^'=^^ if£0(.,t)>O, 

We can understand ( TT3|) as a piecewise heat equation, where the diffusion coefficient 
depends on the convexity or concavity of (p. As pointed out by a referee the equation ( TT3l) 
is studied in [H]. 

We stated Theorem 14.11 for the simple setting of a^ > and / G C^. Theorem 4.6.9 
of [15] states that the Black-Scholes-Barenblatt equation holds for a payoff function with 
linear growth condition: for some a,b > 0, \f{s)\ < a + b\s\, Vs G M. In view of this 
result we expect our Theorem 14.11 also holds for continuous / satisfying a linear growth 
condition. However justifying the limiting argument from discrete time to continuous 
time does not seem to be simple. 



Remark 4.2. The case g^ = needs a special consideration, although Theorem\4-l\ still 



holds for this case. In view of Theorem 4-6.9 of JJEj, the notion of viscosity solution (cf. 
If; IS needed for ^^. 

In Section 6.3 of ',19j this case was treated using parabolic potential theory. The equiv- 
alence of ^3^ to the treatment in Section 6.3 of fWj is seen by the following intuitive 
argument. If q^ = 0, then {d/dt)(j){s,t) > 0, Vs,t, and is non- decreasing in t. strictly 
increases in t at some (so,to) ^/ o-nd only if (f){s,t) is strictly convex in s at this point. 
Then for all t > to, (p{s,t) is (at least weakly) convex in s. This implies that a"^ = a^ if 
and only if(j){s,t) > f{s), which corresponds to the "continuous region" in Section 6.3 of 

The numerical behavior of / [73j) for this case is well illustrated in Figure 6.2 of JT^. It 
should also be noted that (f){s, oo) = limi_>oo 0(s,t) is the least concave majorant (concave 
envelope) of f . 

Note that flT3l) is an additive form of the Black-Scholes-Barenblatt equation in which 
the right-hand side of flT^ multiplied by s^: 

dt 2 ds^ \a^ = q:, if ^0 < 0. 



Consider a multiplicative model, where Reality chooses positive x„'s and Sn = xi x ■ ■ ■ x 
X]y is the product of x„'s. A European option is of the form f{Sj\f). In fl20|) below we 
see that the resulting partial differential equation is exactly the Black-Scholes-Barenblatt 
equation. Note that the multiplicative model is standard in finance literature, although 
it is well known that the pioneering work of Bachelier (j2j|) was formulated in the additive 
form. In this paper we use additive model, because game-theoretic protocols are usually 
formulated in an additive form and also because the limiting partial differential equation 
is a more direct generalization of the heat equation. 

A rigorous proof of our theorem is somewhat tedious and we first give some heuristic 
arguments as to why (TT^ should hold. Later we give a more formal proof, by considering 
an approximate superreplicating strategy as in Section 6.2 of [T9] . 

For our argument it is more convenient to rescale the move space of Market in the 
A^-th game to 

'^N = { I ; • • • , I ; ; — , • • • , , — }• (16) 

After this rescaling, the backward induction in ([8]) for the A^-th game is written as 

/;v(s,Ar-n)=max(pr/^(s + ^,Ar-n-l)+p+J^(s + ^,Ar-n-l)), (17) 

«J VA^ vA^ 



where s = Sn and p^j,p^.i are given by (E]). The initial condition is given by fNiSisf, 0) = 
f{S]\f). Note that by backward induction (ITTj) defines /Af(s, N — n) for all s e M, since / 
is defined on the whole M}. 

As in the proof of Proposition 12.41 fN{s,N — n) is a convex combination of values 
f{s + SN~n)- It should be noted that, unlike the case of convex / in Proposition 12.41 
the weights of the convex combination depend on s. However as seen from the proof of 
Proposition [231 the weights are concentrated either on the two outermost values {a^ , a^} 
or on the two innermost values {aj", a^}, depending on the convexity of /Ar(s, N — n) in s. 
Hence in each interval of convexity or concavity of /Ar(s, A^ — n), it is twice continuously 
differentiable in s. In our numerical studies we found that if the payoff function / is 
smooth and has only finite number of inflection points, then JnIs, N — n) as a function 
of s has no more inflection points than /. 

Write u = N — n — 1. Then each term in the right-hand side of (fT7j) is expanded as 






iv(s + ^,z/) 






fNis, u) + -^/;^(s, u) + R^, (0 < e,, Q, < 1] 
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where derivatives are with respect s and \NRn\ = o(l) uniformly in s and u. Then ((T 
is written as 

Nihis, z/ + 1) - his, u)) = max(-ara+^/;^(s, u) + NRn). (18) 

If we ignore NRj^ = o(l), the right hand side is maximized by {i,j) = {l,Tn) or {i,j) = 
(1, 1) depending on the sign of /j^(s, u). 
Now by rescahng time axis define 

Ms.t) = fNis,Nt), seR,te [0,1]. 

Then ( ITSl) is written as 

This clearly corresponds to (TT3l) . However it seems difficult to let A^ — ;■ oo in (fT9|) and 
prove our theorem directly, although the finite difference approximation to HJB equations 
in Chapter IX of [5] should hold in some form. 

At this point we indicate how the Black-Scholes-Barenblatt equation flTSj) arises in the 
multiplicative case. In the multiplicative multinomial model, Xn is assumed to be of the 
form Xn — 1 G X^, where X^ is given in ( !T6|) . Let S^ = Xi x ■ ■ ■ x x„. Then 

0„_|_1 = 0„ X X„^_i = 0„ + On X [XnJ^x — V). 

The expansion of the right-hand side of ( IT71) in the multiplicative model is 

V-,ifN{s + s\ v) + p+J^(s + s\ v) = Ms, u) + I^s'fUs, y) + Rn- (20) 



This corresponds to ( TT5l) . 

Instead of the above direct approach, knowing that (fT3!) should hold, we can construct 
an approximate superreplicating strategy of Investor and prove our theorem as in Section 
6.2 of [19]. In the following proof, in order to show the inequality 0(0, 1) > limsup^ E{fi\f) 
we adopt a suggestion by a referee. 

Proof of Theorem \4.1\ By Theorem 4.6.9 of ^S\, the solution to (TT^ has a continuous 
first-order derivative in t and a continuous second-order derivative in s. See also Theorem 
11 of [9J and [26]. Consider the following sequence 

0(0, 1), 0(5i, ^^), . . . , 0(5^_i, ^), 0(5^, 0), 

where 4>{Sn, 0) = fiS^)- Writing dSn = Sn+i - Sn = Xn+i, -D„ = 1 - n/N, dD^ = -l/N, 
we can expand the successive difference as 

d(f){Sn, Dn) = (j){Sn+l, -Dn+l) " 0(5'„, -D„) 

d - I d'^ - d - 

= —(t){Sn,Dn)dSn + -^0(S'„, !)„) (rfS'„)^ + —(f){Sn, Dn)dDn + Rn 

= — 0(^„,DOrf'^n - ^(^'(5n,^n) - N{dSnf)—^{S^,D^) + R^, (21) 

11 



where NRn = o(l) is uniformly in s and t. Consider a Markov superhedging strategy V 
(cf. [9]) of Investor which chooses M„ = {d/ds)(j){Sn, Dn). By adding fl2T|) for ra = 1, . . . , A^ 
we have 



g2 



fiSr,) - 0(0, 1) = /C^ - -L ^ (^2(^^^ ^^) _ jv(^^^)2) 7(^^^ ^^) ^ ^(^)^ (22) 



2N^' v"--'M -■v-^-,.. .^^ 

n=l 

At this point we adopt a suggestion by a referee. From Proposition 12.11 and Remark 
12.21 we know that the upper hedging price is computed as the expected value of /{Sn) 
under the extremal risk neutral measure, say p* = p*^. Under any risk neutral measure, 
/C^, n = 1, . . . , A^, is a measure-theoretic martingale and its expected value is zero. Now 
consider the expected value of 

{d^Sn^D^) - N{dSnf)-^4>{Sn,Dn) 

under p* . We can evaluate the expected value, first by conditioning on xi, . . . , x„. By the 
definition of a^, under p* the conditional variance of vNdSn satisfies 

E{N{dSr,Y\x^,...,Xr,)<a^ if ^0(S'„,Z}O>O 

E{N{dSnY I ^1) • • • ) ^n) > cr^ otherwise. 

Therefore taking the unconditional expected value of f l22|) under p* we have E{ffyf) — 
0(0, 1) < except for a term of order o(l). Hence limsup^ E{f^) < 0(0, 1). 

Conversely, consider Market's randomized moves chosen according to the extremal risk 
neutral measure corresponding to E^f^), which is concentrated to two outermost values 
{a'[,a^} or two innermost values {ai,af} at each node of the game tree, depending 
on the sign of 0" (see Remark 12. 2p . Investor's capital is a measure-theoretic martingale 
under this risk neutral measure. Since the measure is supported on two points at each 
node of the game tree, we can modify the standard argument for binomial models to show 
that the expected value of the payoff / converges to 0(0, 1) under the measure. On the 
other hand E{f]\f) is the supremum over all possible moves of Market. Therefore we have 
0(O,l)<liminf^E(/^). D 

Remark 4.3. In the above proof we partly used measure theoretic arguments as suggested 
by a referee. Although we can give a purely game theoretic proof in the line of Section 6.2 
of JTP], it is somewhat tedious. The difficulty lies in the fact that ^(j){Sn, D^) is path- 
dependent. Note that by the game-theoretic law of large numbers (IT^. 01] / ). Investor can 
force that Sj\f/y/N converge to 0. This implies that for large N the empirical distribution 
of Market's moves is approximately a risk neutral measure and Yln=ii^^"-)'^ ^'^ ^^^ variance 
of a risk neutral measure. However because each {dSnf is multiplied by ^(j){Sn, Dn) , the 



convergence Sn/VN — )■ does not imply 



1 ^ r)2 

limmi-J2{^'iSn,Dn)-N{dSnf)j-^iSn,Dn)>0. 
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Although the argument can be fixed by discretization of the values of ^0(5'n,-D„), we 
omit the details. 

Numerically fll3p can be solved by the following backward induction: 1) discretization 
of the interval [0, 1] and M, 2) approximation of the second derivative {d'^/ds'^)(p by the 
second order difference of three neighboring points. Actually this backward induction is 
entirely similar to the exact backward induction in (j8]). When the discretization is not fine 
enough, then the above numerical approximation corresponds to pruning of maximizations 
discussed in Section [31 This suggests that a coarse discretization of the partial differential 
equation yields an approximation which is less than the the true ^(s, 0). 



5 Numerical examples 



In this section we check results of this paper by numerical computation. 

We first calculate the upper hedging price and the lower hedging price of Butterfly 
spread option /(S'at) = max(0, S'at + O.S)— 2max(0, S'at — 0.5)+max(0, Sat — 1.5) inFigure[T] 
under the trinomial model (oi = —l/y/N,a2 = l/vN,a^ = 2/yN). Although Butterfly 
spread does not satisfy the differentiability condition of Theorem 14.11 it can be arbitrarily 
closely approximated by a payoff function satisfying the condition of Theorem 14.11 The 
results are shown in Figure [2] in conjunction with the price under the binomial models. 
From Figure [21 we see that the upper price and the lower price are different from the price 
under the binomial models. 




Figure 1: Butterfly spread 



We now add a new Market's move a^ to the trinomial model and compare the former 
trinomial model (ai = — l,a2 = 1,03 = 2) to this quadnomial model. We consider the 
following three values of 04 as depicted in Figure [31 



1. 04 = -y= (fli < < 02 < 03 < 04) 

2. 04 = ^ (fli < < 02 < 04 < 03) 



"outside". 
"middle". 
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ai --1, a2 = 1, 33 = 2 upper price 

— ai = - 1 , a? = 1 , 33 = 2 lower price 
ai =-1, 32 = 1 

- - a, =-1,33 = 2 




Figure 2: Comparison of binomial model and trinomial model 




^> 




Outside Middle Inside 

Figure 3: Expansion trinomial model into quadnomial model 



q „ _ 0.5 
6. 04 - ^ 



(ai < < 04 < 02 < 03) 



"inside" . 



Figures HI El [6] show the upper hedging prices of the butterfly spread for these quad- 
nomial models compared to those of the trinomial model. In Figure |5] the upper hedging 
prices under the quadnomial model equal those under the trinomial model with increasing 
A^, whereas in Figures HJ [6] the upper hedging prices under the quadnomial models differ 
from those under the trinomial model. 

Next we consider the payoff / = sin(10S'„), which has a lot of changes from convexity 
to concavity, and similarly calculate the upper hedging prices. The results are shown in 
Figures d [H [91 Also in this case the upper hedging prices under the quadnomial model 
equal those under the trinomial model with increasing A^, provided that 04 = 1.5. 

Next, we vary the values of 04 from to 5 by 0.1. Figure [TOI displays the plot of the 
upper hedging prices of the butterfly spread for 1 < A^ < 50 and < 04 < 5. From 
Figure [101 we see that the upper hedging prices converge to an equal value in the interval 
1 < 04 < 2. 
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a, = - 1 , aj = 1 , aj = 2 upper price 
a,=-1,a! = 1, 35 = 2 lowerprice 
ai =-1,32 = 1, 33 = 2,34 = 0-5 upperp 

a, =-1,32 = 1, aa = 2, 34 = 0.5 lower pi 




- 3| = - 1 , 32 = 1 , aa = 2 upper price 
■ a, = - 1 , 32 = 1 , Sj = 2 lower price 

31 = -1, 32 = 1, 33 = 2, aa= 1.5 upper price 
3i=-1, 32 = 1, 33 = 2, 34 = 1.5 lowerprice 










3,=-1, 32=1,3^ = 2 upperprice 

a,=-1,32=1, 3a = 2 lowerprice 

,. ,. a,=-1, 32 = 1, 33 = 2,34 = 2.5 upperprice 
3r=-1, 32 = 1, 33 = 2,34 = 2.5 lowerprice 


Aa^_„.„~^ 




f 



Figure 4: (^ = 04 > 03) Figure 5: (02 < ^ = 04 < 03) 



Figure 6: (-y= 



0.4 < 0,2) 



3i =-1,32 = 1, 33 = 2 Upperprice 
3, =-1.32 = 1, 33 = 2 lowerprice 
3! =-1,32 = 1, 33 = 2,34 = 0.5 upperp 
3, =-1,32 = 1, 33 = 2,34 = 0.5 lowerp 




a,=- 


,32 = 


,33 = 2 upperprice 


--- 3, = - 


.32 = 


,33 = 2 lowerprice 


■ ■ a, =- 


,32 = 


,33 = 2,34=1.5 upperprice 


-- a, =- 


,32 = 


,33 = 2,3, = 1.5 lowerprice 




■ 3i=-1, 32 = 1, 33 = 2 Upperprice 
3r =-1,32 = 1, 33 = 2 lowerprice 
3r =-1,32 = 1, 33 = 2,34 = 2.5 upperprice 
a, =-1, 32 = 1, 33 = 2, 34 = 2.5 lowerprice 




Figure 7: (^ = 04 > 03) Figure 8: (aa < ^ = 04 < 03) 



Figure 9: (-^ = 04 < 02) 



Finally Figure [TT] shows a numerical solution to the partial differential equation ( IT3l) 
for < t < 1 and — 2 < s < 2 for the case of Butterfly spread f{s) = max(0, Sn + 0.5) — 
2max(0, Sn — 0.5) -|- max(0, Sn — 1-5). We compute an approximation of (t){s,t) by the 
following difference scheme: 



m+l 



T^'i A<n 



a 



At 



^i+l 



+ €-1 



As2 



~9 

cr 


= ^^ 


if 0r+i 


-20^ + 0r-i>o, 


~9 


= £^ 


if 0r+i 


-20^ + 0r-i<o. 



(23) 



We rewrite fl23|) as 



in+l 



+ 



a^At,,, 



2As2 



iVi-20r+fc), 



a^ = a^ if 0r+i - 20r + 0r-i > 0, 

0r + 0r-i < 0. 



a^ 



a 



if 



^i+i 



(24) 



We set g^ = I and cr^ = 2 (IT3|) . For discretization we use At = g^^ and As = ^, which 
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price 




Figure 10: Comparison of trinomial model and quadnomial model for various a^ 



satisfies the stability condition (Section 8.4 of [2U] or page 47 of 



At 



1 1 
- < - 
3-2 



for discretization of the heat equation Since our partial differential equation flT3|) can be 
understood as a piecewise heat equation, in our numerical experiments we found that At 
and As satisfying the same stability condition works well. With At = -^ and As = ^, 



we obtain 0(0, 1) ~ 0.3817 and 0(0, 1) ~ 0.2060 using difference scheme f l2l]) . In Figure 
[T2I we compute ( TT4l) for the lower prices by the similar difference method. Table [U shows 
the upper prices and the lower prices obtained in Figure |2l We find that these converge 
to the values obtained by the difference method for the partial differential equations ( 1T3|) 
and (dl]). 



Table 1: Upper prices and lower prices of trinomial model 



A^ 


upper price 


lower prices 


1 


0.2500 


0.0000 


20 


0.3824 


0.1926 


40 


0.3790 


0.1993 


60 


0.3820 


0.2012 


80 


0.3799 


0.2032 


100 


0.3807 


0.2032 
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Figure 11: Numerical solution of the PDE ^ Figure 12: Numerical solution of the PDE dH]) 

6 Concluding remarks 

In this paper we discussed various approximations and asymptotics of upper hedging 
prices in multinomial models. In particular we showed that, as the number of rounds 
goes to infinity, the upper hedging price of a European option converges to the solution 
of an additive form of the Black-Scholes-Barenblatt equation. By numerical experiments 
we checked that this convergence is fast and the asymptotic approximation is useful. 

A multinomial model is the simplest example of incomplete market. A natural exten- 
sion of a multinomial model is the bounded forecasting game ([H]), where Market's move 
is a bounded interval containing the origin. This problem was already considered in [17]. 
Most results of this paper can be extended to the bounded forecasting game. 

Usually the Black-Scholes-Barenblatt equation is studied in the case of vector-valued 
processes. Then the maximum variance and the minimum variance are no longer uniquely 
determined and the maximization in each step of the game tree is more complicated. 
Numerical studies of vector- valued cases are left to our future investigation. 
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